Quenching of pairing gap at finite temperature in ^^^W 
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We extract pairing gap in ^*'*W at finite temperature for the first time from tlie experimental 
level densities of '^^^W, ^^''W, and ^*^W using "thermal" odd-even mass difference. We found the 
quenching of pairing gap near the critical temperature Tc — 0.47 MeV in the BCS calculations. It 
is shown that the monopole pairing model with a deformed Woods-Saxon potential explains the 
reduction of the pairing correlation using the partition function with the number parity projection 
in the static path approximation plus random-phase approximation. 
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Pairing correlations are essential for many-fermion sys- 
tems such as electrons in the superconducting metal, 
nucleons in the nucleus, and quarks in the color super- 
conductivity. The Bardeen-Cooper-Schriffer (BCS) the- 
ory P] of superconductivity has succesfuUy described the 
pairing correlations. This theory was applied to nuclear 
problems at zero temperature ^^r stable nuclei. The 
thermodynamical properties of nuclear pairing were in- 
vestigated by using the BCS theory in the study of hot 
nuclei 0, 0] . Breaking of the Cooper pairs is expected to 
occur at a certain critical temperature in the BCS theory. 

It has recently been reported 1^, ^ that the canonical 
heat capacities extracted from the observed level densi- 
ties in ^'^^Dy, ^'^'^Er and "^Yb form the S shape with 
a peak around the temperature T w 0.5 MeV. These 
S-shaped heat capacities were interpreted as the break- 
ing of nucleon Cooper pairs and the pairing transition 
because this temperature is close to the critical temper- 
ature Tc = 0.57A w 0.5 MeV in the BCS theory For 
the finite Fermi system like a nucleus, however, nucleon 
number fluctuation and statistical fluctuations beyond 
the mean field become large. The fluctuations wash out 
the sharp phase transition, and then the pairing gap A 
does not become quickly zero at the BCS critical tem- 
perature. Several models have taken into account the 
fluctuations beyond the mean field. The quenching of 
pairing correlations have been obtained in recent theo- 
retical approaches: the static path approximation (SPA) 
plus random-phase approximation (RPA) (SL the shell 
model Monte Caho (SMMC) calculations Hi 113, the 
finite-temperature Hartree-Fock-Bogoliubov (HFB) the- 
ory JJj] , and the relativistic mean field theory . 

The odd-even mass difference observed in nuclear 
masses is well known as one of signatures of pairing cor- 
relations. In solid state physics, difference between the 
free energies with odd and even numbers of electrons in 
ultrasmall superconducting grains is found and known as 
even-odd effect 13J. In our previous paper ^^li we sug- 
gested that the suppression of the pairing correlations 
due to finite temperature appears in "thermal" odd-even 
mass difference rather than the S shape of the heat ca- 
pacity. 

In this Rapid Communication, using the thermal odd- 
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FIG. 1: Experimental level densities for ^®^W, ^**W, and 
185-^ The error bars show the statistical uncertainties. 



even mass difference we extract the pairing gap of ^^^W 
from the experimental level densities of "'^^'^W, ^^''W, and 
185^ recently observed [isl It is shown that the 

reduction of the thermal odd-even mass difference is in- 
terpreted as a signature of pairing transition, but not the 
S shape of the heat capacity. 

Figure Q shows the experimental level densities of 
183-w, ^^'^W, and ^^^W extracted from the two-step 7- 
cascade intensities. To study the thermal properties from 
the measured level densities, let us start from the parti- 
tion function in the canonical ensemble with the Laplace 
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transform of the level density p{Ei) 

oo 

Z(r) = ^5i?,p(i?,)e-^'/^, (1) 

1=0 

where are the excitation energies and 5Ei are the 
energy bins. Then any thermodynamical quantities 
0{Z, N, T) can be evaluated by 

oo 

0(Z, N,T) = J2 SE,p{E,)0,e-''^/^/Z{T). (2) 
For instance, the thermal energy is expressed as 

oo 

E{Z, N,T)^Y1 SE,p{E,)E,e-''^/^/Z{T). (3) 

1=0 

The heat capacity are then given by 

We now introduce the thermal odd-even mass differ- 
ence for neutrons defined by the following three-point 
indicator: 

Ai^\z,N,T) = ^-^[E,{Z,N+l,T) 

^2E,{Z, N,T)+Es{Z,N -l,T)l 

(5) 



MeV""'^ and Ei = Ci + A, respectively, where the correc- 
tion factor is given by Ci — — 6.6A~°-^^. The factor / is 
determined so as to adjust the back-shifted level density 
to experimental one. The factors are, respectively, 0.3, 
0.4, and 0.7 for ^^^W , ^^^W, and i^Sw. The parameters 
A for i^-^W and ^*^W are taken as the neutron pairing en- 
ergies 0.61 and 0.72 MeV deduced from mass differences, 
respectively, and for i^**W it is fixed at 1.60 MeV. 

The heat capacities of ^^^W, ^^^W, and ^^^W are 
shown in Fig. 13 All the heat capacities exhibit the S 
shape with peaks around T =0.5 MeV. These heat ca- 
pacities display characteristic behavior similar to those 
of I6i.i62j3y ajjfj i7i,i72Yb 6]. Moreover, we notice that 
the heat capacity of ^^''W around T =0.5 MeV is larger 
than those of ^^^W and ^^^W. The SPA -t- RPA "s] and 
SMMC 01 calculations exhibited the S shape of the heat 
capacity and this odd-even effect where the heat capac- 
ity of an odd-mass nucleus is smaller than that of the 
adjacent even-even nuclei. We can see also the devia- 
tions from the heat capacity of the back-shifted Fermi 
gas model, which is approximated by the Bethe formula 
C = 2aT. In the SMMC calculation O, Liu and Alhas- 
sid identified a signature of the pairing transition in the 
heat capacity that is correlated with the reduction of the 
number of neutron pairs as the temperature increases. In 
their calculation, the pairing correlations are suppressed 
for even-even nuclei, but not for adjacent odd- mass nu- 



where E^ is a shifted thermal energy which is defined by 
substracting the the Coulomb energy from the binding 
energy at zero temperature. The odd-even mass differ- 
ence at zero temperature is known theoretically and ex- 
perimentally as an important quantity to evaluate the 
pairing correlations in a nucleus. The thermal odd-even 
mass difference would be also an indicator of the pairing 
correlations at finite temperature, and is obtained from 
the experimental energies and the level density as well as 
the heat capacity. 

We can calculate the canonical partition function Z{T) 
and the thermodynamical quantities from the measured 
level densities. Formally, the calculations using Eqs. iQJ- 
P|l require infinite summation. However, the experimen- 
tal level densities of Fig. ^ only cover the excitation en- 
ergy up to 6 — 8 MeV. In the evaluation of Eqs. C3l-©j 
therefore, we extrapolate the plots of the experimental 
density to ^ 40 MeV. Here, we use the level densit y fo r- 
mula of the back-shifted Fermi gas model in Refs. [13 
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where the back-shifted energy is U — E — Ei and 
the spin cutoff parameter a is defined through — 
0. 0888 A^/^^/aU. The level density parameter a and the 
back-shifted parameter Ei are defined by a = 0.21A°-^^ 



FIG. 2: Heat capacities as a function of temperature T. The 
broken, dotted, and solid lines denote, respectively, those of 
184-^^ 185^ rpj^g broken dotted line indicates the 
heat capacity of the back-shifted level density formula. 
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FIG. 3: Thermal pairing gap and variation. (a)The thermal 
pairing gap (solid line) extracted from the thermal odd-even 
mass difference Eq. 101 as a function of temperature. (b)The 
variation of the thermal odd-even mass difference defined by 
Eq. (0. 



clei. 

Figure O (a) shows the thermal pairing gap extracted 
from the thermal odd-even mass difference defined by 
Eq. © as a function of temperature. We find a sudden 
decrease of the thermal odd-even mass difference curve 
around T= 0.5 MeV, which is interpreted as a rapid 
breaking of nucleon Cooper pairs and the suppression 
of pairing correlations. We can now regard an inflection 
point of the curve of A^^^ in Fig. |31 (a) as a signature of 
pairing transition, and called it "transition temperature" 
in our previous paper To see more precise position 

of the inflection point, we differentiate A'i\^^'^W) with 
respect to temperature T. It is very important to note 
that the thermal odd-even mass difference has the follow- 
ing identity: 



SAI'^Z, TV, T) 



dT 



{-ir{C{Z,N,T) 



'-[C{Z,N + 1,T) + C{Z,N -l,T)]}. 



(7) 

This identity means that the odd-even difference in the 
heat capacities represents a variation of the pairing cor- 
relations depending on temperature. In Fig. 121(b), we 
can see that the peak of the odd-even difference in the 



heat capacities is a signature of the pairing transition, 
and the thermal odd-even mass difference is a good indi- 
cator for the pairing correlations. The extrapolation of 

f 31 

the level dencity to ~ 40 MeV affects the curve of A„ 
at high temperature in Fig. 3. However, the effects do 
not change the sudden decrease of the pairing gap A„ ' 
around T = 0.5 MeV. 

To describe the above characteristic behavior of the 
heat capacity and the pairing gap extracted from the 
measured level densities of ^^^W, "^*'*W, and "^*^W, we 
consider a monopole pairing Hamiltonian 



H = 



E 

k 



(8) 



where £k are the single-particle energies and P is the pair- 
ing operator P — J2k '^k'^k- By means of the SPA+RPA 
ri9| based on the Hubbard-Stratonovich transforma- 
tion [l8j . the canonical partition function with the num- 
ber parity projection = (1 + cre*''^)/2 (8lll3| IS given 
by 



Tr[P„e 

2 
GT 



^^"^] SPA+RPA 
o 



(9) 
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FIG. 4: Calculated heat capacities and pairing gap as a func- 
tion of temperature T in the SPA-fRPA for the monopole 
pairing model. The upper graph (a) shows the heat capacities 
where the broken, dotted, and solid lines denote, respectively, 
those of ^®^W, ^**''W, and ^^^W. The lower graph (b) shows 
the pairing energy. 
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where a means the even or odd number parity. Here 



•[l + (Tj|tanh2(AfeVT)], 



Crpa = 



a;fcSmh[Afc/^] 
2AfcSinh[cjfc/2r] ' 



(10) 



where Afe = ^if + A^, e'^ ^ £^ _ ^ - G/2, and 7^ = 
£fe ^ M ^ -^fe- The cjfe are the conventional thermal RPA 
energies. If the factor Crpa is neglected, the SPA parti- 
tion function is obtained. Then the thermal energy can 
be calculated from E = —d\nZ^/d{l/T). In this calcu- 
lation, we use the single-particle energies given by an 
axially deformed Woods-Saxon potential with spin-orbit 
interaction |23|. The Woods-Saxon parameters are cho- 
sen so as to fit the experimental single-particle enegies 
extracted from the energy levels of odd nucleus ^'^^Sn 
(^^^Sn core plus one neutron). The deformation takes 
into account effects of a quadrupole-quadrupole interac- 
tion in the mean-field approximation. The deformation 
parameter [3 = 0.23 can be estimated from the exper- 
imental value B{E2) = 119.3 W.u. in the even-even 
nucleus ^^^W. The 50 doubly degenerate single-particle 
energies are taken by assuming ^'^^Sn core, and we fix 
the pairing force strength at G = 20 /A so that the BCS 
pairing gap Abcs reproduce the experimental odd-even 
mass difference A = 0.83 MeV for ^^'^W at zero temper- 
ature. Figure^ (a) shows the calculated heat capacities 
for i^^W, ^^'^W, and ^^^W. We can see that the charac- 
teristic behavior of the extracted heat capacities in Fig. 
121 are well described. Then the neutron pairing energy 
calculated from G{P^P) = GTd\nZ^/dG is also shown 
in Fig. 21 (b). We notice that the calculated heat ca- 
pacity of ^^^W deviates from that of "'^^'^W, in contrast 
with that of -"^^^W, which is different from the result of 



Fig. 2. However, for high temperature this deviation 
would approach to zero according to Eq. (7) because the 

derivative of pairing gap — converges to zero when 
T — > cxD. The monopole pairing model qualitatively ex- 
plains the reduction of the thermal odd-even mass differ- 
ence in Fig. 131 (a) . As pointed out in our previous paper 
[T^ . the peaks of the S-shaped heat capacities is quite 
close to the critical temperature Tc = 0.57A in the BCS 
theory. 

In conclusion, we have extracted the pairing gap of 
^*'*W at finite temperature from the experimental level 
densities of ^*^W, ^^''W, and ^^^W using thermal odd- 
even mass difference. The extracted pairing gap exhibits 
a similar behavior to that of in previous theoretical model 
predictions, that is, the pairing gap decreases from the 
value at zero temperature with increasing temperature. 
The calculations show that while there is no sharp phase 
transition, the pairing gap decreases with increasing tem- 
perature. In particluar, it decreases rapidly around the 
BCS critical temperature. Thus, we can demonstrate 
that the thermal odd-even mass difference is a good in- 
dicator for the pairing transition at finite temperature as 
well as usual one at zero temperature. In this paper, the 
monopole pairing model was used with the SPA-I-RPA 
to describe the heat capacity and the pairing gap, where 
the deformation effect was taken into account. For these 
quantities, however, the fluctuations |3] to the contri- 
butions of the quadrupole-quadrupole interaction should 
be taken into account in more realistic calculations. Fur- 
ther investigations are in progress. We suggest that the 
pairing correlations can be estimated from the measured 
level densities of the triplet nuclei with neutron number 
N + 1,N, and — 1. We hope for further experiments 
to extract the pairing gap at finite temperature. 

One of the authors (K.K.) thanks Dr. A. M. Sukhovoj 
for information of the experimental data and inspiring 
discussion. 
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